COMPUTING 7t(» ANALYTICALLY 

CN ' DAVID J. PLATT 
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^vj ■ Abstract. We describe a rigorous implementation of the Lagarias and Odlyzko 

Analytic Method to evaluate the prime counting function and its use to com- 
r ' pute unconditionally the number of primes less than 10 24 . 

CN ' l. Introduction 
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Computing exact values of the function 7r(x), which counts the number of primes 
less than or equal to x, has exercised mathematicians since antiquity. Early methods 
involved enumerating all the primes less than the target x (using, for example, the 
(— | sieve of Eratosthenes) and then counting them. In 1870 Meissel [T^] described a 

combinatorial method which he eventually used to manually compute tt (lO 9 ) |13j 
(albeit not quite accurately). The algorithm was subsequently improved by Lehmcr 
|11) . then by Lagarias, Miller and Odlyzko [TU] and most recently by Deleglise and 
Rivat [TJ. In 2007 Oliveira e Silva used the algorithm to compute tt (lO 23 ). 
CN , The Prime Number Theorem dictates that all methods reliant on enumerating 

the primes must have time complexity of at least O (xlog -1 x). The latest incar- 
nations of the combinatorial method achieve O [x 2 ^ 3 log -2 x) . 
f~^. ' In their 1987 paper [9], Lagarias and Odlyzko described an analytic algorithm 

l/~) . with (in one form) time complexity O (x 1 / 2+e ). In 2010 Buethe, Franke, Jost and 

£f) \ Kleinjung announced a value for tt (lO 24 ) [5] contingent on the Riemann Hypothesis. 

Their approach "is similar to the one described by Lagarias and Odlyzko, but uses 
the Weil explicit formula instead of complex curve integrals" . This paper describes 
an implementation reverting to Riemann's explicit formula which we have used to 
compute tt (10 24 ) unconditionally. 
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2. A Note on Rigour 

To many, rigorous computation is an oxymoron, due to potential bugs in hard- 
ware, operating systems, compilers and (most likely) user's code. Add to this the 
chance that a power spike or cosmic ray interaction could scupper even a correctly 
written application, and the situation seems hopeless. We do what we reasonably 
can to minimise such risks including running applications on systems with ECC 
memory after testing on hardware from different vendors with different operating 
systems and using different compilers. 

However, there are certain aspects over which we do have more control. Esti- 
mating the rounding error that will accumulate through a complex floating point 
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computation is a non-trivial task that we eschew. Instead we rely on interval arith- 
metic (see jH] for a good introduction). Thus, instead of storing a floating point 
approximation to any real value, we hold an interval that brackets the true value 
with the standard operators and functions being overloaded to handle such inter- 
vals. 

Furthermore, we will estimate some quantities by, for example, truncating an 
infinite sum. We will need to derive a rigorous bound for the error introduced, but 
rather than keep track of such errors manually, we can simply add them to the 
interval being evaluated and let the software take over. Indeed, we can in some 
circumstances use interval arithmetic to compute such bounds for us. 

3. The Analytic Algorithm 

The analytic algorithm relies, at its core, on Perron's formula. 

Theorem 3.1 (Perron's Formula). Let a(n) be an arithmetic function with Dirich- 
let series 

a(n) 



fl(-) = E 



n" 

n=l 

Now if g(s) is absolutely convergent whenever 3?s > a a , then for c > a a and x > 
we have 

c+zoo 



E* fl («) = ^- / g(s) 

*■ — ' Zm J 



ds 
x" — 
s 



where the * on the summation sign indicates that if x is an integer then only 1/2 
of the a(x) term is included. 

Proof. See page 245 of [T] and the subsequent note. □ 

The relevance of Perron's formula to the matter at hand comes from the series, 
absolutely convergent for 3fts > 1, 

lft , n f A(n) 

log COO = 2^ 



n s log n 

n=2 b 



Here A is the von-Mangoldt function so -r-^- takes the value — at prime powers 

° log n m L L 

p m and zero elsewhere. Define 



tt*(x):= J2 ~ 



m 

p rn <x 



where if x is a prime power we only take 1/2 of its contribution to the sum. Then 
applying Perron's formula we get for c > 1 and x > 

e+'ioo 

(3.i) £^ = .*(..) = -L /bgcOO^. 

*■ — ' logn Zm J s 

n ^ x c-ioo 

Here we note that, although we can cheaply recover ir(x) from 7r*(:r), the slow rate 
of convergence of the integral dooms any attempt to use it in this context. 
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At this point, Lagarias and Odlyzko introduce a "suitable" Mellin transform pair 
4>{t) and 4>(s) and derive 

(T + 'ioo 

(3-2) **(x) = ± I \ogas)M s )d S + J2-lx x (p m )-^(p m )} 

a— ioc ^ 

where Xx{t) is the characteristic function of x, defined by 

f 1 t<x 

Xx(t):=\ 1/2 t = x. 

{ t> x 

We note that taking 4>(s) — — makes cf)(t) = Xx{t) and we recover Equation l3.il 
Thus estimating tt*(x) now splits into estimating an integral and summing 4>(t) 

evaluated at prime powers in the vicinity of x. 

In his PhD thesis [8] Galway investigated the proposed algorithm and suggested 

using the Mellin transform pair 

-. . x s /A 2 s 2> 

4>{s) := — exp 

and 

0(t) := 



V2X 
Here erfc is the complementary error function 




erfc(x) := — = / exp (-t 2 ) At 

V7T J 



and A is a positive real parameter used to balance the convergence of the integral 
with the width of the prime sieve. 

Galway showed that <f> and <f> as defined are indeed "suitable" and, using argu- 
ments based on the uncertainty principle, suggested that they are in some sense at 
least optimal. He also gave a rigorous bound for the error introduced by truncating 
the prime sieve to some finite width. 

a+ioo ^ 

4. Evaluating ^7 / logC(s)<A(s)ds 

cr — zoo 

At this point, we depart from the line taken by Lagarias and Odlyzko. Rather 
than attempt to numerically estimate the integral in Equation 13. 2\ we take an 
approach closer to the spirit of Riemann and evaluate it in terms of the non-trivial 
zeros of (, leading to Theorem 14.71 below. 

Before proceeding, we need a couple of lemmas. 

Lemma 4.1. The "Round the Pole" lemma. Let f be a meromorphic function with 
a simple pole at a with residue R, and let T be the semicircular contour anticlockwise 
from a + e to a — e. Then 

lim / f(z)dz — iriR. 

e^0+ J 

r 
Proof. See page 29 of [H]. □ 
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Lemma 4.2. We have 



Proof. 



lim(logC(l + e)-log(-C(l-e)))=0. 

e->0+ 



lim(logC(l + e)-log(-C(l-e))) 



limlog^A 



= lim log 



l/e + 0(l) 



o+ l/e + 0(l) 



D 



Lemma 4.3. There exists a sequence ofTj — > oo such that for any a G [— 1,2] we 
have for s = a + iTj 

C -{s) = 0{\og 2 T 3 ). 

Proof. Referring to Davenport [6], for any zero /3 + 17 of C with 7 large, there 
are 0(log7) zeros with imaginary part G [7—1,7+1] (Corollary (a), page 99). 
Therefore we can select a Tj within 0(1) of 7 such that Tj differs from the imaginary 



part of any zero by 3> 



logTj 



By (4) on page 99 we have for a G [—1,2] 



jfr + iT/, 



E 



a + iTn 



■OQogTj 



where the sum is taken over zeros with imaginary part G [Tj — 1,7} + 1]. There 
are C(logTj) such zeros, each one making a contribution to the sum limited by 
©(logTj) and the result follows. □ 

Lemma 4.4. For t G K 

|log(-C(-l + *t))|<5 + t 2 . 
Proof. By the functional equation 

C(-i + ^ = c(2-,t)i-ggI^I^). 



We then use 



' 1+it ) T ( 1+rt ) = r (^-5^) so we can apply Stirling's approximation. 
Also, for 5Rs > 1 we have 



|IogC(*)l = 



E 



A(n) 1 



log(n) n s 



< 



E 



A(n) 



1 



^ log(n) n 



SF =logC(Ks). 
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Lemma 4.5. 



-1+ioo 



— / log(-C(*))^(*)da 



exp(f) / 2 



Proof. We use Lemma 14.41 and take absolute values, majorising with 



exp 



(*) 



2ttx 



5 / exp ( — — - dt + / |i| exp ( — - — ) dt 



where both integrals can be evaluated. 



□ 



Lemma 4.6. Let $(s) be the unique holomorphic function $ : C \ R<o — > C such 
that 



• <J?'(s) = 0(s) and 

$(<r + it) + $(a - it) 



• lim 



= for any fixed real a . 



Then 



(1) 3?(s) — logs extends to an entire function, 

(2) lim $((7 + it) = C is purely imaginary and 

t— >oo 

(3) <&(er ± it) =F C is rapidly decreasing as t — > oo. 
Proof. To show (1) we define for s ^ K<0 

s 

F(s) = / l(z)dz 
l 
where the contour of integration is the straight line from 1 to s. We then define 

F(l + iT) + F(l-iT)~ 



$(s) := lim 

T— >oo 



F(s) 



and we have 



F(s) — logs 



CO-- |<1:. 



To show (2) we observe that $(s) = $(s) and from the definition we have C+C = 0. 
To show (3), we take T > and we have 



\$(a±iT) T C\< ^ eX p(^f- 



exp 



2.,.2 



-A 2 t 



dt 



D 



Theorem 4.7. Let <&(s) 6e as defined in Lemma \476\ Then 

2+zoo — 1+ioo 

X - J 0(s)logC(s)ds = $(l)-^5R$(p)-log(2) + ^ J 0(s)log(-C(s))ds. 



2?ri 



Proof. We will refer to the contours represented in Figure [TJ 
These contours are 
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2+ioo 



Figure 1. Contours to evaluate -^ J 4>(s) log C(s)ds 



2 — ioo 



r 2 
r 3 
r 4 
r 5 
r 6 
r 7 
r 8 
r 9 



the semi-circle clockwise from 1 — e to 1 + e for e small and positive. 

the semi-circle clockwise from 1 + e to 1 — e. 

the horizontal line from 1 + e to 2. 

the horizontal line from 2 to 1 + e. 

the vertical line from 2 to 2 + iTj, Tj not the ordinate of a zero of (. 

the vertical line from 2 — «T 9 to 2. 



the horizontal line from 2 



iTj to -1 



iTj. 



the horizontal line from — 1 — iTj to 2 — iTj. 



the vertical line from —1 



iTj to 



followed by the clockwise 



circular arc centred at —1 to 
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the clockwise circular arc centred at — 1 from i to — 1 — 



followed 



by the vertical line to 



-1-iTi 



the horizontal line from i to 1 — e. 
Ti2 - the horizontal line from 1 — e to j. 
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We consider the integrals 

(4-1) ^ /($(*)- C)^ds 

2m J Q{s) 

for the contours I\, T3, Ts, Tj Tg and Tn in the upper half plane and 
(4.2) JL /($(,) + C )^d. 

for F2, T4, T6, Tg Tio and Ti2 in the lower half. 

We denote the integrals in Equations 14.11 or 14.21 as appropriate along T n by I n 
and proceed as follows. 

For J 5 and I 6 we get 



lim (J 5 + I e ) = lim r— : 

J— >oo ?'— >oo Z7TI 



= lim : 

i— >-oo 27ri 

1 
~ 2?ri 



$(«) - C) ^( s ) ds + J (*( s ) + C) ^(s)ds 
Lr 6 r 6 . 

*(*) - C) logC(s)[ +i '' + (*(») + C) logC(s)' 2 



2-iT; 



(s)logC(s)ds 



1 
2tt7 



2+ioo 

2ClogC(2)- /" ^( S )logC(s)d S 

2 — ioo 

Considering the contours T 7 and r 8 , we use Lemma |431 and the Gaussian decay 
of <&(s) ± C from Lemma T4. 61 to conclude 



lim (J 7 + h) = lim — - 

j— >oo j— >oo Z7TZ 

= 0. 



$(s) - c) ^(s)ds + / ($(s) + c) ^(s)ds 
lr 7 r 8 ' 



Considering Ig and 7io we have 
1 



lim (ig + Jio) = lim . 

J-»oo J — S-oo ^7Tl 



= lim - — r 

j— S-oo 27TI 



*(*) - C) ^(s)ds + /" (*(«) + c) ^(s)ds 

-To Tio 

1/4 



*(*)-C log(-C(a)) 



($(*) + c) log(-C(*)) 



-1— tT. 

1/4 



1 

2?r7 



1 

"2vr7 



(«)log(-C(*))da+ / 0(a)log(-C(*))dfl 

Lr 9 r" 



(a)log(-C(*))d* + 2Clog(-C(l/4)) 



ToTio 
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where the contour of integration is Tg followed by Tio- Convergence of this integral 
is due to Lemma l4~5l and the zero free region of £(s) with \s + 1| < | and $ls > — 1. 
For In and I12 we have 



hi +It2 = 



2ni 

1 

2^1 



S(s) - c) j{s)ds + f ($(a) + C) ^(s)ds 

T11 T12 

($(«) - c) log(-C(s))r £ + ($(*) + C) log(-C(s)) 



1 

2?ri 



■(a) log(-C(s))ds + / 0(a) log(-C(s))ds 

iTn ri 



1/4 



l-£ 



— [2Clog (-( (1/4)) - 2Clog(-C(l - e))] 



For h and I2 we find 



h+h 



2ni 

1 

2^1 



LTi 



(*(*) - C) ^( S )d, + / (*(«) + C) ^( S )d 



C, 



r 2 



Vi r 2 



27TI 



| (s)ds _y|( s)ds 

LTi r 2 



by Cauchy's Theorem since the residue of 4- at 1 is — 1. 
Finally, for 7 3 and Z4 we get 



h+h = 



2ni 

1 

2?ri 



($(.s) - C) C -{s)ds + / ($(«) + C) ^( S )ds 



lr 3 



*(«) - C ) log C(s) + ( *(«) + C ) log C(s) 



l+e 

2 



1 

2~7ri 



(s)logC(s)ds+ / 0(s)logC(s)ds 
lr 3 r, 



— [2ClogC(l + £ )-2ClogC(2)] 

27TI 



Now by Cauchy's Theorem and exploiting the fact that the non-trivial zeros of 
C, occur in complex conjugate pairs, lim J2k=i h = 12 ^^(p) so we have 



COMPUTING tt(x) ANALYTICALLY 



E«*w = $ « - h 



2+'ioo 



a)logC(a)ds+ / <f>(s) log(-{(s))ds 



2— too 



r B ,ri 



%ogC(l + e)-log(-C(l-e))] 



C 

2^i 



Of) 

coo 



ds 



Of) 

coo 



ds 



Now the result follows from taking the limit as e — > + by Lemmas 14. 21 and 14. II and 
then straightening the line of integration of the second integral to Sfts = — 1 . This 
introduces a contribution of log(— £(0)) = —log 2 from the pole of 4>(s) at s = 
with residue 1. □ 

Again, if we take <f>(s) = — then $(s) = Ei(logs) where Ei is the exponential 
integral and we recover Riemann's explicit formula 

-l+zoo 



7r*(a;)=Ei(logx)-^Ei(ploga ; )-log2+ f log(-C(s))x f 



ds 

s 



We truncate the sum over zeros so we need a rigorous bound for the error that 
this introduces. We derive such a bound in Appendix [A"l 
The computation of n(x) now reduces to 

• enumerating the prime powers near x, 

• computing <f>(t) at these prime powers, 

• locating the non-trivial zeros of C to sufficient accuracy and 

• evaluating $ at these zeros (and at 1). 

5. The Prime Sieve and cj>{p) 

To compute ir (10 24 ) with the zeros at our disposal we needed a sieve of width 
rj 6 x 10 15 centred at 10 24 . We will discuss only locating the primes in that interval, 
the prime powers being a trivial task by comparison. 

Two basic methods were considered, sieving (necessarily segmented) and a hy- 
brid technique described by Galway [S] . The latter proceeds by first eliminating all 
y-smooth numbers and then applying a base 2 Fermat primality test to the remain- 
der. Given a list of the (few) numbers in our range which are composite, y-rough 
and yet still pass the Fermat test, we are done. Our tests suggest that while an 
implementation of the Hybrid sieve would not be competitive at height 10 , the 
crossover might not be far away. 

Our implementation used Atkin and Bernstein's sieve based on binary quadratic 
forms [2] to enumerate the sieving primes (< x 1 ' 2 ) feeding a segmented implemen- 
tation of the sieve of Eratosthenes to delete composites in the target region. 
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y t 1 

y = a 3 x 3 / 

/ El 

aw / J 


Y' 


X 



Figure 2. Approximating a Cubic with a Line (Lemma 15.1 



For each sieve segment centred at xq, we output 

•EL 

p 

• J2( x o ~P) and 
p 

• E<>o-p) 2 - 

By restricting the sizes of the segments, we can ensure that the entire computation 
can be achieved using native 64 bit integer instructions, with the exception of the 
third sum which requires 128 bit addition. However, this represents only a small 
performance penalty on modern CPUs. These three terms are then used to form an 
approximation to ^2 0(p) by Taylor series. In fact, three terms are not enough to 

p 
give us the required precision so we exploit the following lemma to derive a linear 

approximation to the fourth (cubic) term as well. 

Lemma 5.1. If we approximate the real cubic y — a 3 x 3 on the interval x € [—w, w] 

o 2 

where w > with the line y = ax with a — 3 ™ , then the magnitude of the error 



over the interval is < 



Ml 



What is more, in terms of minimising the worst case 



error, this line is the best choice of any quadratic. 

Proof. We refer to Figure [5] Without loss of generality, take 03 > and since both 
a^x 3 and ax are odd, we consider only the interval x € [0,w]. The error E\ is 
simply 03 w 3 — aw and E2 is at its maximum where the slopes of the line and the 
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cubic are equal. This happens at x = ^ hr— so _E? = \/-S- \ / tM — . The worst 

^ ^^ y 3a3 V 3a3 V 7as 

case error follows from setting E\ = E2 and solving for a. 

The maximum error occurs 4 times at x G {±w, ±, /g 5 -}- This means that any 

curve which improves on the line must be below the line at x £ {— w, . g^ - } and 

above it at x € |w, —./■^ t -\. Thus, such a curve would have to cross the line at 
least 3 times which is not possible for a quadratic. □ 

6. Computing $ 
The following lemmas give us a means of computing JR$(/o). 
Lemma 6.1. For 3fe 7^ and feeR 

ex p( ~ A 2 2/l2 ) 
0(s o + «/i) = <f>(so) cxp (i/i(s A 2 + log(a;))) 7~TJh ■ 

so 

Proof. We start with 

cxp / A 2 (s +i/t) 2 ^ ^sa+ih 
4>(s + lh) = — 

s + in 
and rearrange to get 

exp (^) x S0 cxp (ih(s X 2 + log(a;))) cxp l - x '~"~ 



so 1 + ^ 

so 



Lemma 6.2. Let N E 2Z >0; A, h > and Ah < 1. Then 

\ / „ =0 V / 

where 

1 /A 2 h^ f 



a 



1^1 < 



(f)! 



Proof. This function is entire and the restriction on Aft, makes the terms alternating 
in sign and decreasing. □ 

Lemma 6.3. Let N <G Z>o and \h\ < \so\- Then 

ih\ \—\ ( —ih\ „ 

i + - =J2( — + Eb > 



where if R = 



have 

R N 



\Eb\< 



1-R 



Proof. This function is analytic on the open disk \h\ < |so| and the missing terms 
form a geometric series. □ 
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We can now fix some N <G 2Z>o and multiply these two (degree N) polynomi- 
als to yield a single (degree 2N) polynomial in h which we can integrate against 
exp(ih(X 2 + log(x))) analytically. 

We now start at <& (i) and move up the | line in small steps. We take the 
real part of the highest p used to be zero and bound the error this approximation 
introduces. 

7. The Computation 

To obtain the O (a; 1 / 2 " 1 " 6 ) time complexity of Lagarias and Odlyzko's algorithm, 
we should choose the free parameter A to equate the run times of the sum over 
zeros and the sieving elements of the computation. In fact, we biased the run time 
towards computing zeros, both to confirm RH holds to a height sufficient for this 
computation and since this data may have application to future researchu Wc 
isolated all the zeros of ( to a height of 30, 610, 046, 000 (103, 800, 788, 359 zeros). 
The technique used to locate these zeros is described in [T5] but is in essence a £ 
specific, windowed version of Booker's algorithm from [4]. 

We set A = 6273445730170391 x 2~ 84 (note that this is exactly representable in 
IEEE 754 double precision floating point). We used the first 69, 778, 732, 700 zeros 
to compute the sum (those to height 20, 950, 046, 000) which in turn dictated that 
we sieve a region of width about 6 x 10 15 . As a consequence, the truncation error 
from summing over the zeros and from the sieve were together < 0.989. 



With this choice of A, we have 



*(§) 



< 3 x 10 13 so we need our zeros to be 



located to an absolute accuracy of at least 25 decimal placesi 2 ] Thus, we are forced 
to use multiple precision arithmetic, despite the performance penalty this implies 
(up to a factor of 100 compared with hardware floating point). 

As discussed earlier, we use interval arithmetic to manage the accumulation of 
rounding errors during the computation and to this end we have extended Revol 
and Rouillier's MPFI package [16] in the obvious way to handle complex arithmetic. 
Adopting interval arithmetic incurs another performance penalty (a factor of about 
3 or 4 this time). 

The sieving (entirely in integer arithmetic) was performed on 352, 800 segments 
each of width 2 34 as dictated by memory constraints, further sub-divided to control 
the error in approximating <fi by Taylor series. The actual computation of this 
Taylor approximation again requires multiple precision interval arithmetic. 

The sum over zeros and the prime sieve all parallelise trivially and we used the 
University of Bristol Bluecrystal Phase II cluster to perform all the computations, 
consuming approximately 63, 000 CPU hours. In a personal communication Tomas 
Oliveira e Silva indicated that computing 7r (lO 23 ) using the combinatorial method 
required about a month on a single computer, so it appears that we are still some 
way short of the crossover implied by the asymptotic run times of the two algo- 
rithms. 



To this end, Jonathan Bober has made about the first 36 billion zeros available at [3]. 
Our zeros are located to an absolute accuracy of ±2 _lt)2 which more than suffices. 
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The result of the computation, after adding in the various error terms, was an 
interval straddling a single integer, so we have 



Theorem 7.1. 



7T (10 



24\ 



18, 435, 599, 767, 349, 200, 867, 866. 



We note that this agrees with the conditional result of Buethe, Franke, Jost and 
Kleinjung. 



Appendix A. Truncating the Sum over Zeros 

We require a rigorous bound for the error introduced by truncating the sum over 
zeros in Theorem 14. 71 We proceed as follows. 

For t > 0, t not the imaginary part of a zero of £, define N(i) to be the number 
of zeros of £(s) in the critical strip with 3s £ [0, i\. 



Lemma A.l. Lett>2. The 



n ^-tM^ L 



< 0. 137 log(i) +0.443 log log(i) + 1.58 



Proof. See Q2] Theorem 19. 

Lemma A. 2. For x > 1, T, A > and a € [0, 1] define 



D 



B(a,T) :=exp 



A 2 (l-T 2 



1 



Tioga; X 2 T 2 x 



Then 



5R$(cr + iT) <B(a,T). 



Proof. We integrate along the contour running vertically down from — 1 + ioo to 
— 1 + iT, then right to a + iT. For the horizontal contour we have 



exp 



( \ 2 (u+iT) 2 \ 

\ 2 / ^u+iT 



u + iT 



For the vertical contour we have 



du 



< 



"-(yo /,. 



exp 



< 



(ggg) 

Tioga; 



/ A 3 (-l+<t)' \ 

-1 + zt 



^ exp (^i±^ 



,-1+it 



df 



< x 1 exp 



exp 



-x-'r 



( \ 2 \ oo 

exp(-j , 
< -^ / '-P 



-A 2 £ 2 



di 



dt 



exp 



/ X 2 (l-T 2 ) \ 



A 2 T 2 



D 
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Lemma A.3. Let T > 0, a £ [0, 1] and a T be such t aT > N(t) for all t > T. The 



Y^ B{a,^p) <exp 



•3p>T 



A 2 (l-T 2 ) 



1 



Tlogx X 2 T 2 x 



Proof. Referring to Lemma IA.2I and writing k a '■= exp I — ) 
can majorise the sum with the Stieltjes integral 



[ A2T2 + 2 -^V(T)l 


■) 


x" | 1 


T\ogx ' X 2 T' 2 x 



wo 



\ 2 t 2 " 
£; CT exp( — - — ) <UV(t). 



We now integrate by parts and majorise N(t) with t aT to obtain 

oo 

-A 2 T 2 \ ,„ k„ 

y(c, ^p) < — k a exp I - 

3p>T ^ 



£ B(a, Qfp) < -K exp ( -^ j JV(T) ^ ( 



X t A exp 



-A 2 i 2 



d* 



— k n 



X 2 T 2 + 2 

\2j^2-a T 



exp 



-X 2 T< 



exp 



-A 2 T 5 



N(T) 



a 



We note that the olt referred to above can be computed using Lemma [A. II 
We now consider the error introduced by truncating our sum over the zeros of 

C- Let T\ be height below which we find and use all the zeros, and Ti the height to 

which we know RH to hold. 

Lemma A. 4. Let E\ be real part of the error introduced by ignoring zeros with 
imaginary part of absolute value £ [Ti, T2] (whose real parts are all known to be ^)- 
Then 



I £1 1 < 2 exp 



A 2 (l-T 2 



1 



Tilogx X 2 T^x 



A 2 2\ 2 + 2 
X 2 T^° Tl 



N(Ti) 



Proof. We apply Lemma IA.3I with a — | and introduce a factor of 2 for the zeros 
with negative imaginary part. □ 



We note this bound includes all the zeros with imaginary part > T2 but their 
contribution will be negligible. 

Lemma A. 5. Let E2 be the real part of the error introduced by omitting the zeros 
with imaginary part $ [— T%,T2\ from the main sum, (which do not necessarily have 
real part = -k ). Then 



I-E2I < exp 



A 2 (l-T s 



x + 1 2 

Talogx + X 2 T 2 x 



X 2 T 2 + 2 



a 2 t; 



^--^2) 



Proof. We pair each p with 1 — p and take the worst case which is when one of the 
zeros has real part 1 and the other 0. The result then follows from Lemma [A. 31 □ 
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